Breakdown of the operator product expansion in the 't Hooft model 
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We consider deep inelastic scattering in the 't Hooft model. Being solvable, this model allows us to 
directly compute the moments associated to the cross section at next-to-leading order in the 1/Q 2 
expansion. We perform the same computation using the operator product expansion. We find that 
all the terms match in both computations except for one in the hadronic side, which is proportional 
to a non-local operator. The basics of the result suggest that a similar phenomenon may occur in 
four dimensions in the large N c limit. 
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It has been more than 35 years since QCD was 
vindicated as the theory of the strong interactions 
1] . Unfortunately, analytic solutions describing the 
hadrons and their interactions in terms of the de- 
grees of freedom (quarks and gluons) and parameters 
(coupling constant and quark masses) that appear in 
the QCD Lagrangian have remained elusive. Leav- 
ing aside symmetry considerations (or how they are 
realized), the only quantitative and analytic scheme 
to check the dynamics of QCD from first principles 
is through (perturbative) weak-coupling computa- 
tions. The problem then is to relate those computa- 
tions to experiments, which are naturally described 
in terms of hadrons. In order to do so two key in- 
gredients are used: 

1) Asymptotic freedom [l|, Q, which states that 
the renormalized coupling goes to zero at large Eu- 
clidean momentum. 

2) The operator product expansion (OPE) Q and, 
more generally, factorization. So far, the OPE has 
only been proven within perturbation theory [J]. 

In practice, the combination of these two points is 
stated as the fact that the multiplication of two op- 
erators (sandwiched between physical states) enjoys 
the following expansion at short distances 



A(x)B(Q) = ^C n (x)O n (0) 



(1) 



where O n are local operators with increasing dimen- 
sionality in n and the coefficients C n can be com- 
puted in perturbation theory (or at least at weak 
coupling). 

The (Fourier transform of the) OPE probes these 
correlators at large Euclidean momentum. There- 
fore, the OPE is not directly accessible to experi- 
ment and one has to resort to dispersion relations 1 
to check it. This practical version of the OPE is 



at the basis of computations at large Euclidean mo- 
mentum of (the moments in) deep inelastic scatter- 
ing (DIS) and the vacuum polarization tensor, which 
so far have been thought to be among the more solid 
predictions of QCD, since they are not affected by 
quark-hadron duality problems 2 . Therefore, the im- 
portance of setting the OPE and the factorization 
methods used in quantum field theories, specially 
in QCD, on solid theoretical ground can hardly be 
overemphasized. So far, it was thought that the use 
of the OPE (in its non-perturbative formulation [fjj) 
was secure, even tough it has not been proven in 
QCD. It has been only partially checked in models, 
for instance in two dimensional QCD in the large 
N c limit: the 't Hooft model pj. This theory is su- 
perrenormalizable and asymptotically free, so it is a 
nice ground on which to test the OPE 3 . This was 
done at the lowest order in the OPE in Refs. @, Q 
for the vacuum polarization and for DIS off a me- 
son with nice agreement between the results of the 
model and the OPE expectations. In Ref. [1J| the 



the lack of knowledge of the asymptotic behavior of the 
correlators to ensure that one can neglect the contributions 
at infinity when using the Cauchy theorem. This will not 
affect the main conclusions of this paper, since this ambi- 
guity should have an impact on, at most, a finite number 
of moments. 

We do not enter in this paper on the use of the OPE and 
factorization methods for quantities living in, or affected by, 
the Minkowski cut. These quantities are usually regarded 
as less fundamental, and are the ground on which the dis- 
cussion about quark-hadron duality takes place, see for in- 
stance [5| and references therein. For those observables one 
can easily find examples where perturbation theory fails 
in the large N c limit. For instance, the imaginary part of 
the vacuum polarization tensor becomes a sum of infinitely 
narrow resonances in the large N c limit, as opposed to the 
smooth result obtained from perturbation theory. 
In the 't Hooft model there are no marginal operators. 
Therefore, the coupling constant does not run and has di- 
mensions; no renormalons should then arise. 



1 These may have potential problems of their own; basically 



OPE was numerically checked at next-to-leading or- 
der (NLO) in the 1/Q 2 expansion, with logarithmic 
accuracy, for the vacuum polarization. In this paper 
we consider DIS at NLO. We will only give the main 
results and postpone the details to a later publica- 
tion. In practice, we will compute the correlator 

T^(q)=i [d 2 xe^(z3;n\Tf(x)f(0)\tj;n) (2) 



pM 






/'/: -^K^)t(q 2 ,x b ) 



where xb = Q 2 /(2P n ■ q), Q 2 = —q 2 , P n is the mo- 
mentum of the meson and i, j stand for the fla- 
vor of the quark and antiquark, respectively, which 
form the bound state. j 7 '(a;) = J2h3h( x )> wnere 
Jh( x ) = e h 4> h ^il) h (x). 

The imaginary part of T^ v is proportional to the 
differential cross section when xb > 0, and in the 
light-cone frame reads 



ImT++ = 



^ 



dP+ 



2^7 2(2tt)P,.. 

x(2n) 2 5 2 {q + P n -P l 



(ij;m\j + {0)\ij;n)\ 



"E 

<) I Ml 



(ij;m\ ^2 e h jt(0)\ij;n) 

h=i,j 

- 2 (1 -x B ) 



Mi 



■I'B 
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Using dispersions relations, up to a finite number of 
subtractions, one should have (see also the discus- 
sion in Ref. [ill ]) 



T{Q 2 ,x B ) = - 



where 



1 ImT(Q 2 ,y B ) ( .. 
d VB~ ; — — , (4) 
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x = 1/(1 + (P 2 - P 2 )/Q 2 ). Therefore, 
T(Q 2 , xb) admits an analytic expansion in 1/xb for 
x B > x^ ax , 



T(Q 2 ,x B )=A J2 M N (Q 2 )-^ , 



(5) 



Af=0,2,4, 



where 



M N (Q 2 ) ee 1- 
Zir 



dy B y%- l ImT(Q 2 ,y B ). (6) 



General expressions for the matrix elements in terms 
of the 't Hooft wave functions, 4^{x), will be pre- 
sented in Ref. 12]. They are obtained using sim- 
ilar techniques to those in Ref. [lj], where the 
Lagrangian, the 't Hooft equation, as well as the 



definitions of the bound states are given. Approx- 
imated expressions, valid when 1 — xb ^> 1 IQ 2 
(f3 2 = g 2 N c /(2n) is the 't Hooft coupling), will 
also be obtained in Ref. [12| using Eqs. (51) and 
(52) from [lCj], as well as similar techniques to those 
used in Ref. |R|. For the explicit computation it 
is convenient to work with the kinematical variable 
X = —q + /P+ (we take q + < and q~ — ► +oo), 

which satisfies the equality xb — x/(l — -j^-x 2 ). 
One then obtains ImT with 0(1/ Q 2 ) precision when 
l-x>/3 2 /Q 2 
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"i n. — m 1~ P 2 - Up to a prefactor, the terms 
ei/jTUi/j {■ ■ ■} represent the contribution from the 
matrix element of the current of the quark i and 
the antiquark j, respectively. Note the relative sign 
— (— l) m between both contributions, which can be 
obtained from symmetry arguments. 

Eq. ([7|) is one of the key results of our paper. By 
inserting this expression in the moments and using 
the Euler-MacLaurin expansion we obtain (x max = 
1 - M 2 /Q 2 + 0(1/Q 4 )) 
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which is correct with 0(1/Q 2 ) precision (for finite 
N). Note that with this precision we can replace 
Xmax = 1- Note as well that the last two terms in Eq. 
dHJ) are 0(m 2 /Q 2 ) suppressed with respect to the 
leading term. One comment is in order here. When 
using the Euler-MacLaurin expansion we had to deal 
with terms proportional to (—1)'™. Such terms give 
a potential contribution of relative order o(l/Q 2 ), 
beyond the precision of our computation. Neverthe- 
less, we also have terms proportional to (— l) 2m = 1. 
They come from the interference of the particle and 
antiparticle currents and are no longer suppressed by 
the sign-alternating behavior, only by the prefactor, 
which is 0(l/Q 2 ). Therefore, these terms contribute 
at 0(1/Q 2 ). We anticipate that they will be the ones 
responsible for the violation of the OPE. 

Using Eq. {3} and the approximated expressions 
we have obtained for the moments in Eq. ©, we can 
obtain an approximated expression for T, which can 
actually be written as a dispersion relation formula: 



^(^)- 2 r 1 4y B toT»*-{<P,VB) 
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where ImT EuL is given by Eq. (|7|), performing the 
substitution ^ m — ► J dm. T EuL is not a good ap- 
proximation to T when we approach the physical 
cut. This has to do with the fact our expressions for 
the moments are not valid when N — > oo (x — > 1 
limit). Nevertheless, it can be considered a generat- 
ing functional for the moments with not very large 
N. Therefore, it is the expression we would expect 
to be equal to the OPE result. In order to ease the 
comparison, we write T in a factorized way: 



T EuL (Q 2 ,x B ) 



Q 2 



(10) 



x / dy{e?J i (x,y)f i (y)+e 2 j J j (x,y)f j {y) 
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In the above expressions we have not inserted the 
Wilson line: 



*(aT,y-) = P[e (i9f y- dz A+{z )} ] 



(14) 



between the quark fields ip to make gauge invariancc 
explicit, since we are working in the light-cone gauge, 
A+ = 0. 

The functions J are defined as 4 
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(16) 



We can now perform the same computation us- 
ing the OPE (we need to do the computation at one 
loop). In order to avoid spurious differences with 
the hadronic result, we compute ImT OPE and use 
dispersion relations afterwards (actually, one can al- 
ready see the failure of the OPE calculation from 
the comparison of ImT OPE and ImT Eul ). We then 
obtain 

T OPE (Q 2 ,x B ) = -2^ 2 (17) 

dy{e 2 J i (x,y)f i (y)+e 2 J J (x,y)f j (y)} . 



4 It is possible to redefine J so that it has the functionality 
J(y/x). This implies redefining f(y) by some powers of y, 
equivalent to introducing some extra 9+ derivatives. We do 
not do so in this paper because it would increase the length 
of the formulae. Note as well that the derivatives in J; tend 
to change the variable of fi(y) to fi(y(l + m^ R /Q 2 )), con- 
sistent with the interpretation from the perturbative com- 
putation in the OPE. 



We can easily see that T OPE ^ T EuL . If we consider 
the moments generated by T OPE , M® PE , they can 
be expressed in terms of non-perturbative local op- 
erators (for the explicit expressions see pjl) as ex- 
pected. Nevertheless, Mn has extra terms. If we 
consider the difference, we obtain 



M N (Q 2 ) - M° PE {Q 2 ) = 16 ^i - "':"'< 



Q 4 N c Q 2 



(-iD ) N + 2 
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x^(z-)^l + (0)^ + (z-)\i 3 :n) , (18) 

which is expressed in terms of non-local operators. 
We take this result as evidence of the existence of 
OPE-violating terms. Let us stress that this is the 
first time, that we are aware of, that an analytic 
calculation in a quantum field theory exhibits OPE- 
violating contributions. 

The possible existence of OPE-breaking effects in 
QCD has already been discussed in the past. As 
early as in Ref. [lj] numerical evidence for the exis- 
tence of OPE-breaking effects in the gluon conden- 
sate was claimed. Nevertheless, it is still unclear 
whether those effects can be associated to ultravi- 
olet renormalons and/or higher orders in perturba- 
tion theory (for a recent discussion see [la]). Over 
the years there has also been some discussion on the 
possible existence of a (A 2 ) min condensate. This ob- 
ject should actually correspond to a non-local gauge- 
invariant condensate, though its explicit form is un- 
known for QCD [16] . Finally, there are some models 
that may produce effects that break the OPE, see 
for instance [13|. Nevertheless, those OPE-breaking 
effects would affect the static potential and the vac- 
uum polarization. Regarding this we would like to 
emphasize that we do not find any OPE-breaking 
effect in the static potential or the vacuum polariza- 
tion in the 't Hooft model. The static potential can 
be computed exactly in the 't Hooft model within 
perturbation theory Therefore, there is no room 



there for effects associated to a sort of (A 2 ) m i n , con- 
densate. With the present precision of our compu- 
tation, we also do not sec OPE-breaking effects in 
the vacuum polarization. 

The existence of this OPE-breaking effect is much 
associated to the large N c analysis we have done, and 
the existence of sign-alternating effects. It should 
not be difficult to devise similar large iV c -inspired 
models in 4 dimensions, which would produce OPE- 
violating terms. In those models, however, it might 
be difficult to disentangle the OPE-violating terms 
from standard OPE contributions, since they both 
would scale in the same (or a very similar) way and, 
at the end of the day, one fits higher twist effects to 
data. 

The new term is analytic in 1/Q 2 and can actu- 
ally be written in a factorized form. Therefore, one 
may think whether one could extend the standard 
formulation of the OPE in order to include this sort 
of terms. It would also be interesting to get a bet- 
ter understanding of this term from a diagrammatic 
analysis, if possible. 

In conclusion, we have shown that in the 't Hooft 
model, at NLO in the 1/Q 2 expansion, the moments 
associated to DIS receive a contribution that: 

1) cannot be written in terms of local operators, 

2) cannot be matched with any OPE-like contribu- 
tion we are aware of. 

We take these two facts as "smoking-gun" signals 
for the breakdown of the OPE in DIS in the 't Hooft 
model at NLO in the 1/Q 2 expansion. Note that it 
appears as a subleading (NLO) effect in DIS. More- 
over, it is important that we considered the inter- 
ference between two currents, otherwise this effect 
would be suppressed by a factor of, at least, order 
1/Q 4 . For the vacuum polarization we do not see 
this sort of effects with the present precision [lCJ • 
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